Destructed double-layer and ionic charge separation near the oil-water interface 
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We study suspensions of hydrophobic charged colloidal spheres dispersed in a demixed oil-water 
mixture by means of a modified Poisson-Boltzmann theory, taking into account image charge effects 
and partitioning of the monovalent ions. We find that the ion's aversion for oil can destroy the 
double layers of the oil-dispersed colloids. This affects the salt-concentration dependence of the 
colloidal adsorption to the oil-water interface qualitatively. The theory also predicts a narrow range 
of the oil-dielectric constant in which micron-sized water-in-oil droplets acquire enough charge to 
crystallize at volume fractions as small as ~ 10~^ in the absence of colloids. These findings explain 
recent observations [M.E. Leunissen et al., Proc. Nat. Ac. Sci 104, 2585 (2007)]. 



PACS numbers: 82.70.Kj, 89.75.Fb, 68. 05. -a 

Colloidal particles can strongly adsorb to an oil-water 
or an air-water interface, and hence form a monolayer. 
Since a pioneering study by Pieranski jj] a lot of atten- 
tion has been devoted to the lateral structure and the 
lateral colloidal interactions of such monolayers 0, 
This two-dimensional picture is often realistic because 
the "binding" potential of a single colloid to the in- 
terface can be very strong, typically of the order of 
10'^ — lO^ksT for a micron-sized colloid due to surface- 
tension or image-charge eflFects [l|, 0, Here ks is 
the Boltzmann constant and T the temperature. How- 
ever, in very recent experiments by Leunissen et al. [^, 
on systems of charged, micron-sized, hydrophobic poly- 
methylmethacrylate (PMMA) spheres dispersed in an oil- 
like mixture of cyclohexyl bromide and cis-decalin, the fo- 
cus was not only on the in-plane structure of the adsorbed 
colloidal monolayer at the oil-water interface, but also 
on the out-of-plane "perpendicular" colloidal structure. 
Striking observations of this three-dimensional study in- 
clude (i) a dramatic increase of the lattice spacing of the 
oil-dispersed colloidal crystal up to 40 /im by bringing 
this dispersion in contact with water, (ii) strong colloidal 
adsorption to planar and spherical oil-water interfaces 
even for non-wetting colloids, (iii) an extremely large 
colloid-free zone between the adsorbed monolayer and 
the bulk crystal in the oil phase, in some cases with a 
thickness d > 100/im, and (iv) the existence of micron- 
sized water-in-oil droplets that are sufRciently charged 
to crystallize into a lattice with a spacing of 10-15 mi- 
crons, without any colloid presence in the system, pro- 
vided the dielectric constant of the oil was in the nar- 
row range 4 < Co ^ 10 0|. In this Letter we develop 
a theory for a mixture of hydrophobic colloids, cations, 
and anions in an oil-water mixture, taking into account 
screening, image charges, and self energies. We do not 
focus on the lateral structure [J, [3, 0] but instead on 
the structure perpendicular to the interface, and not on 
a single colloid [s, 0| but on a many-body system. Key 
new mechanisms that we identify this way are (i) double- 
layer destruction, due the preference of ions to reside in 
the water phase rather than in diffuse layers surrounding 



the oil-dispersed colloids, and (ii) ionic charge separa- 
tion if the cation's and anion's aversion for oil differs, 
such that water-in-oil droplets can acquire a net charge. 
These phenomena explain essentially all observations of 
Ref.0|, and could be ingredients to further understand 
and manipulate Pickering emulsions 0, S, S|- 

We consider a planar water-oil interface in the plane 
X = 0, characterized by a macroscopic surface tension 
-fwo- The interface separates two semi-infinite continu- 
ous bulk phases of water {x < 0) and oil {x > 0). This 
system is a medium for a three-component mixture of hy- 
drophobic colloids (radius a, charge Ze), cations (radius 
a+, charge e), and anions (radius a_, charge — e). Here 
e is the elementary charge. The strongly hydrophobic 
character of the colloids is described here phenomeno- 
logically through the colloid-water and colloid-oil surface 
tensions jcw = lOmN/m and jco = ImN/m, respectively. 
Note that these tensions do not incorporate electrostatic 
contributions. Following Pieranski's geometric argument 
[l]] a colloidal particle with its center at a; G (—a, a) is 
therefore subject to the external potential 

X x^ 

V{x) = 27ra^(7c^ - 7co)(l ) - 7ra^7^o(l r), (1) 

a 

while a particle completely immersed in water (x < —a) 
or oil {x > a) has V{x) = 4:na'^{-/cw - Jco) — lO'^/csT 
for a ~ 1/im and V{x) — 0, respectively. Note that we 
shifted the potential of Ref.[l| by an arbitrary constant 
for later convenience. The potential V{x) has a deep 
minimum at a; = a;* = a{'jcw — Jco)/jwo = a cos 9 with 9 
the wetting angle, provided \ jcw ~7co\ < "two- Otherwise 
V{x) is monotonic and we speak of non-wetting. Below 
we consider the wetting case 7^,0 = 9.2mN/m, such that 
cos 6* = 0.987 and V{x*) ~ -lO^fc^T, from which strong 
adsorption of micron-sized colloids at x ~ is expected. 
We also consider the non- wetting case "iwo = 9mN/m 
such that cos 9 — 1, where no strong adsorption is to be 
expected because V{x) is monotonic. 

The planar water-oil interface also generates an exter- 
nal potential for the ions due to the dielectric discontinu- 
ity, which leads to different electrostatic self energies in 
the two solvents. We write the self energy in medium i — 
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w,o as the Born energy e^/(2eia±) = f±{ei). Denoting 
the local dielectric constant by the step function e{x) = 
Cw for X < and Cq for a; > 0, this self-energy effect can 
be accounted for in terms of (conveniently shifted) ionic 
external potentials V±{x) = /±(e(x)) — f±{ew), which 
vanishes in water and is of the order of (1 — 20)A;bT for 
realistic eo — 4 — 20 in the oil phase, i.e. the ions prefer 
to be in water. 

Now that we have specified the external potentials 
V{x) and V±{x) for the colloids and the monovalent ions, 
respectively, we employ the framework of density func- 
tional theory to calculate the equihbrium density profiles 
p{x) and p±{x) [ill, [l3, [l3|. The grand-potential 
functional p+,p_] is written, per unit lateral area 
A, as 

j= JdxpUx)[kBninP^ - l) + VUx)) 

+ J dxp{x) (fcBT(ln ^ - 1) + V{x)) 

+ ksT J dx[p{x)^{fjix)) + ^Qix)cPix)), (2) 

where ri{x) = Ana^ p{x)/3 is the colloidal packing frac- 
tion, and where the first and second line are the ideal- 
gas grand-potential functionals of the ions and the col- 
loids in their external fields, respectively, and the third 
line describes the hard-core and Coulomb interactions. 
The chemical potentials are represented in terms of an 
ion concentration ps and a reference colloid packing 
fraction 770 to be discussed below. The colloid-colloid 
hard-core interactions are taken into account by the 
Carnahan-Starling excess free energy per particle 5" (77) = 
(477 — 377^)/(l — 77)^ with the weighted packing fraction 
f]{x) = J^^2a dx'w{x — x')r]{x') with the (low-density) 
weight function w(a;) = -^{Aa^ — x'^)/a^ [l3|. Such a non- 
local treatment of the hard-core interactions is necessary 
to describe the extremely localized adsorbed colloidal 
monolayer in the case of wetting realistically. The elec- 
trostatic interactions between all species are described 
in Eq. |[2]) at a mean-field level in terms of the total local 
charge number density Q (a;) — Zp{x)+p+{x)~'p^{x) and 
the yet unknown electrostatic potential kBT(j){x)/e that 
must satisfy the Poisson equation and boundary condi- 
tions 

e(x)(/)"(x) = -AnPe^Qix), (x 7^ 0); (3) 
lim (x) = limeo^'(x) ; lim (/)'(x) = 0. 

x'iO xlO x^±oo 

Here a prime denotes a derivative with respect to x, and 
(3 = l/{kBT). Note that the second line of ^ de- 
scribes the image-charge effects (first term) and global 
charge neutrality (second term). Minimizing the func- 
tional leads to the Euler-Lagrange equation 6Q./5p{x) = 



50./ 5p±{x) = 0, which can be rewritten as 

77(2;) = 7/oexp(-/3y(x) -Z0(x) -/3/i(x)); 
p±{x) = p,cyix>{ - I3V±{x)t<P{x)), (4) 

with I3p{x) = ^{n{x)) + dx'w{x-x')r]{x')^' inix')) 

where ^''(77) = d'i' {rj) / drj . Note that the strongly hy- 
drophobic character of the colloids leads to a vanishing 
colloid-density at x ^ —00, such that 0(— 00) — and 
P±{~oo) — Ps for electro neutrality reasons, i.e. the bulk 
water phase acts as the salt reservoir with a total ion 
concentration 2ps. The corresponding Debye screening 
length in the water phase is = {SnPe^ ps / ew)~^^^ ■ 
In the bulk oil-suspension, x co, the average colloid 
packing fraction 7y(cxD) = 77^ can be imposed by tun- 
ing 770 appropriately. Eqs.|[3l) and ^ also imply, for 
X > 0, that (j)"{x) = Klsmh{(j){x)-(t)c)-4:nl3e'^Zp{x)/eo, 
with the Debye screening length in the oil phase de- 
fined by kI = K^(e.u,/eo) exp[-/3(V+(oo)+l/_(cx)))/2] and 
<^, = /3(F_(oo)-y+(oo))/2. 

Eqs.llHl) ^nd lUl) form a closed set for the four unknown 
profiles p{x), p±{x), and <p{x), and can be solved nu- 
merically on an x-grid by standard iterative methods on 
desktop PC's. Typically we use several hundred non- 
equidistant grid points, with a relatively small grid spac- 
ing close to X = and close to the monolayer position 
X = X* in the wetting case. 

We fix as many parameters as possible in accordance 
with those in the experiments described in Ref.[^, such 
that Z = 450, a = 1pm, r]b = 5x 10"^, e„ = 80, = 5.2. 
For simplicity we take equal ionic sizes a+ — — 0.3 
nm such that V±{x) = ITfcsT for x > for both ionic 
species. This leads to /a = 1 — 40 if one varies ps 
from 0.07 to 26 mM. In Fig. [l]we show, for these param- 
eter choices, the resulting colloidal packing fractions 77(x) 
for several screening lengths k~^, in (a) for the nonwet- 
ting case cos 6 = 1 and in (b) for the partial wetting case 
cos6 = 0.987. In all cases the colloids are so hydrophobic 
that 77(x) is vanishingly small for x < 0. The qualitative 
difference between the two sets of curves in (a) and (b) 
is striking, given the small difference in contact angle. 
In Fig. 1(a) we observe 77 (x) to increase by up to several 
orders of magnitude near the interface at a < x < 5a, 
whereas in Fig. 1(b) we observe a densely packed mono- 
layer of colloids at X = X* separated from the bulk by a 
colloid-free zone of a thickness d of the order of several 
K~^. Although the increase of 77(x) up to 10^776 close to 
the interface in Fig. 1(a) is significant, it is much weaker 
than to be expected on the basis of the attractive image 
charge potential W{x) = Z Z'e'^ /{Acqx) ~ — 500afcsT/x 
that a single colloid (without cations and anions) would 
experience at x > a for the present parameters, with 
Z' = —Z{ew — f-o)l{(-w + £0) = —395 the image charge 
[5]. At first sight one would attribute this relatively mod- 
est colloidal adsorption to the mutual repulsion between 
adsorbed colloids, which in the present many-body de- 
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Figure 1: The packing fraction profile 77(2;) of strongly hy- 
drophobic, oil-dispersed colloidal spheres (radius a = 1/im, 
charge Z = 450) in the vicinity of a planar interface at 
a; = between water {x < 0, dielectric constant = 80) 
and oil {x > 0, to = 5.2), for a colloidal bulk packing frac- 
tion rji, — r;(oo) = 5 x 10~^, in (a) for nonwetting col- 
loids {cos 9 — 1) for oil screening lengths K~^/a = 4 — 40 
from bottom to top, and in (b) for weakly wetting colloids 
(cose = 0.987) for K^'^/a = 4-40 from left to right. The 
insets show, for = 40a, the salt and colloidal charge distri- 
butions, revealing the complete (a) and partial (b) destruction 
of the double layer around the adsorbed colloidal layer. 



scription is included through the mean-field Coulombic 
term in Eq.(l2]). However, in that case one would expect 
a stronger (longer-ranged) colloid-colloid repulsion and 
hence a weaker colloidal adsorption for increasing 
whereas Fig. 1(a) shows the opposite trend. This implies 
that another mechanism is at work here. The mechanism 
that we identify involves the salt and colloidal charge dis- 
tributions Qsix) = p+(x) — P-{x) and Qc{x) = Zp{x), 



respectively, shown in the insets of Fig.l for K~^/a = 40. 
These insets show a charge separation due to the accumu- 
lation of a net negative ionic charge in water, in a very 
thin layer — < x < 0, and a net positive colloidal 
charge for a < x < 5a in (a) and at a; = a;* in (b). In 
other words, the usual double-layer structure, in which 
the compensating ionic charge resides in the vicinity of 
several Debye lengths of the colloidal charged surfaces, is 
destroyed completely (a) or partially (b) for hydrophobic 
colloids in the vicinity of the oil-water interface: a frac- 
tion of the compensating ionic charge prefers to reside in 
the water phase "further away" rather than close to the 
colloidal surface. Since the length scale for local charge 
separation is set by the Debye length, this fraction, and 
hence the net charge in the thin layer — < x < in 
the water, will increase with increasing k~^. This accu- 
mulated charge attracts the oppositely charged colloids 
at a; > a in Fig. 1(a), and hence the adsorption of colloids 
increases with k~^. In a an oversimplified dynamic sense 
one could envisage a colloid with a spherical double layer 
approaching the water-oil interface from the oil side, then 
being stripped from (part of) its ionic cloud if it gets to 
close, such that it cannot diffuse back to the bulk as a 
neutral entity and hence adsorbs. Double-layer destruc- 
tion also takes place in the wetting case of Fig. 1(b), where 
the deep well of V{x) at x — x* drives the formation of a 
close-packed monolayer aX x — x* . Due to colloid-colloid 
repulsions this monolayer strongly repels the colloids in 
bulk, giving rise to the colloid-depleted zone of a thick- 
ness d that increases with expected on the basis 
of this argument, in agreement with Fig. 1(b). However, 
whereas d ~ 4k~^ for < 10a, d is levelling off for 
larger to a maximum value of about d ~ 80a at 
= 50a, beyond which d is decreasing with (not 
shown). This pecuHar behavior of d at large is en- 
tirely due to the double-layer destruction, which causes 
(an increasing fraction of) the compensating charge of 
the monolayer at a; = a;* to reside in the water phase 
at — < X < 0, such that the geometry tends to be 
that of two oppositely charged plates which have a strong 
field inside but a weaker field outside, i.e. the repulsion 
between the monolayer and the bulk colloids is reduced 
and hence d shrinks. This non-monotonic behavior of d 
with is a mere prediction of the present theory, but 
the observations of d decreasing from about 100 to 20 mi- 
cron upon adding salt in the water phase (3| (and hence 
varying over the parameter regime of Fig. 1(b)) is in 
agreement with the theoretical values. 

One problem in the comparison between the results of 
Fig.l and the experiments of Ref.[^ is that the colloid- 
depleted zone was observed for non- wetting colloids, 
whereas Fig. 1(b) holds for the case of wetting colloids, 
cos 61 = 0.987. One could on the one hand of course ar- 
gue that the wetting angle of the experimental system is 
smaller than the resolution of its measurement, such that 
the experimental system would actually be (extremely 



4 



weakly) wetting. But on the other hand we found in ex- 
tended calculations that profiles such as those of Fig. 1(b) 
can also be obtained for non-wetting colloids provided 
we take into account the combined effects of (i) charge 
regulation, such that an additional profile Z{x) is calcu- 
lated using ionic association-dissociation equilibrium as 
an additional condition, and (ii) a nonlocal analogue of 
the local relation Qc{x) = Zp{x), such that a colloid at 
X contributes to the charge density in the whole interval 
[x — a, x+a\. These extension will be pubHshed elsewhere. 

In Ref.^] it was also reported that water-in-oil droplets 
could be stable without any colloids or any other addi- 
tives in the system. This observation goes against the 
common believe that emulsions require "emulsifiers" in 
order to be stable The stabilization mechanism 

proposed in Ref.[l| is based on the asymmetry between 
the cations and anions as regards their self-energy in 
oil and water, such that the droplets spontaneously ac- 
quire a net charge. Moreover, it was observed that a 
system of water-in-oil droplets could actually even crys- 
tallize, but only if 4 < < 10. The present theory as 
formulated in Eqs.(l3])-(l4l) can be employed to underpin 
and further understand these recent surprising observa- 
tions. The absence of any colloids is modelled by setting 
p{x) = throughout, and the alleged asymmetry of the 
cations and anions is taken into account by setting their 
size ratio a_|_/a_ = ^ < 1, such that V-{x) = S^V+{x), 
i.e. the cations' aversion for oil is larger than the an- 
ions'. For ^ = 0.5 and e„ = 80, we solved Eqs.®-!]!]) 
on an x-grid to obtain the equilibrium profiles p± (x) , for 
Ps[-^] = 10~^, 10~^, for a range of eg and a+. The net 
(positive) charge in the water phase (per unit area, in 
units of e) follows as cr = J_^dx{p+{x) — p_(x)). If 
we now presume that a spherical water-in-oil droplet of 
radius R has the same surface charge density, it would 
have a total charge number Z — AirR^a, and two of these 
droplets would repel each other by a potential of the form 
vir) = (Ze)^ eyii,{Ko{2R - r))/[eoil + KoR)^r] In- 
spired by the experiments of Ref.ffi] we set R — 1.5^m, 
and consider Pv{r) for the inter-droplet spacing r = 10 R 
as a function of Eq and a+ as shown in the contour plot 
of Fig. 2. It is clearly observed that v{10R) only exceeds 
IksT in a very narrow band of values for eo. Given that 
one expects crystallization of a dilute system of charged 
spheres as soon as v{r) at typical droplet spacings exceeds 
O^ksT), the narrow eo-bands of Fig. 2 correlate with the 
observations in Ref.@|. The peak in Z at Co — 30 ^ 10, 
shown in the inset of Fig. 2, shows that v{10R) is deter- 
mined by an intricate balance of amplitude and range; the 
fact that ZXb/R < 5 for eo < 10, with As = (3e^/to the 
Bjerrum length of the oil, indicates that charge renormal- 
ization is of no concern in the regime where l3v{\QR) ~ 1 

In conclusion, we have presented a theory for charged 
colloids, anions, and cations near a water-oil interface, 
taking into account colloidal wetting properties, ionic 



self-energies, electrostatic image effects and screening. 
We focussed on the structure perpendicular to the in- 
terface, and identified the destruction of colloidal dou- 
ble layers and the separation of ionic charge as im- 
portant mechanism taking place near this interface. 
With these non-lateral mechanisms we can describe sev- 
eral recent observations of hydrophobic charged colloids, 
most notably the large colloid-depleted zone in between 
the adsorbed monolayer and the bulk oil-dispersed col- 
loidal phase, the spontaneous self-charging of water-in-oil 
droplets, and their crystallization in a narrow band of 
Q. An extension of the theory to describe the very re- 
cent observations of stable Pickering emulsions [l5| is in 
progress. 
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fully reading the manuscript. 




10 20 30 

^0 

Figure 2: Contours of the pair potential v{r) between water- 
in-oil droplets of radius R — 1.5/im at separation r = 107?, 
as a function of cation radius a+ and dielectric constant of 
the oil Eo, for two water salt concentrations ps ~ lO/^M and 
O.IM. The anion radius is O- = 2a+ here, and no colloidal 
particles are present. Within the similar narrow region of eo 
we find that v{im) - v{9R) = C»(fcsr).The inset shows, for 
a+ — 0.3nm and ps = 10~^M, the charge Z and combination 
ZXb/R of the droplet (see text). 
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